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Abstract. Recently it has been shown that when the Dark Matter (DM) particles captured 
in the Sun directly annihilate into neutrino pairs, the oscillatory terms in the oscillation 
probability do not average to zero and can lead to a seasonal variation as the distance 
between the Sun and Earth changes in time. In this paper, we explore this feature as a novel 
method to extract information on the properties of dark matter. We show that by studying 
the variation of the flux over a few months, it would in principle be possible to derive the DM 
mass as well as new information on the flavor structure of the DM annihilation modes. In 
addition to analytic analysis, we present the results of our numerical calculations that take 
into account scattering and regeneration of neutrinos traversing the Sun. 
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1 Introduction 

If the dark matter is composed of Weakly Interacting Massive Particles (WIMPs), it can be 
trapped inside the Sun and give rise to a neutrino flux which is in principle detectable at 
the neutrino telescopes. Detection of such a flux is an indirect way of detecting DM which 
has received much attention in the recent years. If the flux is high enough, in addition to 
establishing WIMP as the DM, information on DM annihilation can be derived from the 
properties of the flux such as energy spectrum or flavor composition. 

Recently it has been shown in [1] that if the annihilation to neutrino pairs is one of the 
dominant modes, the oscillatory terms in the oscillation probability will lead to a seasonal 
variation of the number of muons produced by the Charged Current (CC) interaction of 
Ufj^ from the Sun in the neutrino telescopes. That is because in this case, the spectrum of 
neutrinos is monochromatic so, contrary to the general belief, the oscillatory terms do not 
average to zero and vary throughout a year as the distance between the Sun and Earth changes 
due to nonzero eccentricity of the Earth's orbit. As pointed out in [1], the seasonal variation 
can be considered as an independent tool to derive information on the DM annihilation 
modes. In this paper, we further elaborate on this possibility. We analyze the information 
that can be derived from this observable on the flavor structure of DM annihilation amplitude. 
We also present numerical results using a code that solves evolution equations describing the 
neutrino propagation, taking into account the effects of neutrino absorption, scattering and 
fr regeneration inside the sun. 

Our particular attention is given to the case that DM particles primarily annihilate into 
neutrinos. As has been discussed in detail in ref. [2], wide classes of models can be built in 
which the dominant annihilation mode of DM particles is the annihilation into neutrinos. At 
first sight it might seem that within such models, the DM interactions with nuclei will be too 
weak to give rise to a substantial DM abundance in the Sun center and thus to a significant 
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neutrino flux. We briefly discuss this issue and show that, despite dominantly annihilating 
into neutrinos, DM can have large enough spin-dependent scattering cross section off the 
protons inside the Sun. Our main approach in this paper is however model independent. Our 
goal is to find out to what extent the properties of DM (e.g., Br(DM + DM — )■ faJ^is)) can be 
determined by combining the information on seasonal variation and on flavor composition 
(more precisely, the ratio of the detected muon-like events to shower-like events). 

This paper is organized as follows. In section 2, the general properties of the flux is 
described and the seasonal variation is quantified. In section 3, the theoretical framework 
within which a DM pair dominantly annihilates to a neutrino pair is discussed and shown 
that the DM capture rate in the Sun and the subsequent neutrino flux can be high enough to 
make the method presented in this paper viable. In section 4, the properties and limitation of 
neutrino detectors are described and the observable quantities are formulated. In section 5, 
the numerical code that has been developed to carry out the calculation is described. In 
section 6, the information that can be in principle derived from the observable quantities 
defined in the previous sections on the flavor structure of DM interactions are analyzed. In 
section 7, numerical results are presented and the observed patterns are analyzed. Concluding 
remarks are given in section 8. 

2 Neutrino flux from dark matter annihilation in the Sun 

DM particles propagating in the space between stars and planets in the galaxy can enter the 
compact objects such as the Sun or the Earth. If the scattering cross section of DM off nuclei 
is large enough, these particles can lose energy at scattering and fall into the gravitational 
potential of the Sun or the Earth. As a result, the DM particles will be accumulated in the 
center. The accumulated DM particles annihilate with each other and produce the Standard 
Model (SM) particles. Among the SM particles only neutrinos can reach the surface from 
the Sun or Earth center. The rest will either be trapped or decay before reaching the surface. 
DM annihilation can give rise to a neutrino flux either directly (i.e., DM + DM — t- I'u) or 
as secondaries (i.e., DM -|- DM — t- XX and subsequently X — t- I'Y). If the mass of the DM 
particles is above a few hundred GeV, the neutrino flux from the DM annihilation can be 
detectable at the neutrino telescopes such as IceCube. The number of the neutrino events at 
IceCube produced by DM annihilation in the Sun center can be of the order of a few hundred 
events per year. The background neutrinos pointing towards the Sun is below 100 events 
per year [3-5] , so the detection of an excess of these neutrinos at these energy ranges can be 
interpreted conclusive indirect detection of DM. 

Depending on the couplings and characteristics of the DM particles, the flavor com- 
position and the energy spectrum of the neutrino flux can be different. In particular, if 
neutrinos are secondaries, we expect a continuous spectrum. On the contrary, if DM par- 
ticles annihilate directly into neutrinos, the neutrino energy spectrum at production point 
will be monochromatic. That is because the DM particles annihilate non-relativistically. In 
fact, the thermal motion of DM particles can widen the line to a very narrow Gaussian with 
width 5E/E ~ 10~^ [1]. Propagation of the monochromatic neutrinos from the dark matter 
annihilation in the Sun has been studied in a number of papers including in [7-9]. This 
difference in spectrum can in principle be invoked to discriminate between the scenarios pre- 
dicting different decay modes. Ref. [6] has systematically studied the possibility to extract 
Br(DM + DM — t- z^i/) and Br(DM + DM — ?• rf ) by measuring the energy spectrum. 
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The spectrum of neutrinos emerging from the Sun surface will consist of a sharp line, 
corresponding to un-scattered neutrinos, superimposed on a continuous tail corresponding to 
scattered or regenerated neutrinos. The sharp line can still give rise to oscillatory behavior 
which is the base of the method suggested in this paper to determine if (DM + DM — )• z^z/) is 
among the dominant annihilation modes. Following ref. [1], we define the seasonal variation, 
A, as 

Afa.Ati;t..At.).^""^'''-^"''^'-) (2.1) 
JV(ti,Ati)+iV(t2,At2) 

where 

N{t, At) = — ^ 2.2 

j;-'^'A,ff{e[t])L-^t)dt 

in which L is the distance between the Sun and the Earth and Aeff(6[t]) is the effective area 
of the detector which depends on the angle between the beam direction and the axis of the 
detectors, 0(t). As the Earth orbits around the Sun, L and 6 both change with time during 
a year. Notice that with this definition, the dependence of the flux does not affect A. 
A vanishes when the oscillatory terms are absent or averaged out. As shown in [1], A can be 
substantially large and exceed 50 % so with a moderate accuracy, its nonzero value, which 
is a proof of the existence of the mode (DM + DM — )• vu) can be established. The value 
of A depends on the DM mass, Br(DM + DM — )• vu) and the initial flavor composition as 
well as neutrino oscillation parameters. As pointed out in ref. [1], the value A can be used 
to extract information on the initial flavor composition, especially when it is combined with 
energy spectrum measurements. Of course, to construct A, the chosen intervals, Ati and 
At2, should not be too small. Otherwise, the statistics will be too low, making the derivation 
of A meaningless due to the large statistical error. For the energy interval of our interest 
{i.e., El, ~ 100 — 500 GeV) even if we take Ati + At2 = 1 year, the oscillatory terms do not 
average out. Even the oscillatory terms corresponding to Am|^ survive. We will elaborate 
on this point later on. Considering this fact, there is practically no upper bound on Ati, 
except that Ati + At2 < 1 year. 



3 Theoretical motivation 



As discussed in ref. [1], in order to have non- vanishing A, the following two conditions are 
necessary: (1) The neutrino spectrum has to include a sharp line with width 5E/E <^ 
A-kE / {Am? L) which can be realized if the DM annihilation into neutrino pair is significant; 
(2) the neutrino flux at production must be non-democratic; i.e., F^^ : F^^ : F^^ ^1:1:1. 
We have already discussed the reason for the first condition. The second condition is also 
necessary because for F^ = F^^ = F^^ = F^^ , the neutrino flux at Earth will be independent 

of the oscillation probability {F^^ = Y.a^uc.^i'^a ^ ^Z?) = ^ ^^z?) = -^°) and 

as a result, A will vanish because it comes from the oscillatory terms in the oscillation 
probability. 

As discussed in ref. [2], various classes of models can be built that satisfy both these 
conditions. An explicit example is given in ref. [10]. In order to measure A, the number 
of collected neutrino events has to be larger than a few hundred. If the flux is close to the 
upper bound from AMANDA, such an amount of data can be collected at IceCube. For 
n^DM ~ 100 GeV, the neutrino emission from the Sun with a rate of order of 10^° sec~^ 
will be enough to lead to a few hundred events per year. In the saturation limit where the 
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capture rate of DM particles by the Sun, Cq, equals the neutrino emission rate (i.e., two 
times the annihilation rate), the neutrino emission rate can be translated into a bound on 
the DM-nuclei cross section. As seen from eqs. (9.25) and (9.26) of [11], for spin-independent 
scattering mediated by a scalar, with a DM-nucleon cross section, cr„, as small as 10~^ pb, 
Cq = 10^*^ sec~^ can be achieved. Such cross section is below the present bound on spin- 
independent DM-nucleon cross section [12]. It is noteworthy that while the Sun is mainly 
composed of protons, for spin-independent interactions, the heavier nuclei inside the Sun play 
the dominant role in capturing the DM particles. That is because the nucleons in a nucleus 
interact coherently with DM so the cross section grows quadratically with nucleus mass. If 
the DM-quark interaction is mediated via a heavy scalar, we would have 



Thus, (cj(DM + DM — <?g)v) [decoupling ~ 10"^ ph <C {(Jtotv) ~ 1 ph. The above discussion 
shows that it is possible to construct a model within which DM-nucleus interaction is large 
enough for sufficient DM capture in the Sun but at the same time, the main annihilation 
mode of DM pairs is the annihilation into neutrino pairs. 

4 The IceCube neutrino telescope 
4.1 General description 

IceCube is a km'^-scale neutrino telescope recently completed at the south pole [13]. IceCube, 
which encompasses the AMANDA experiment [14], has been designed to detect the flux of 
galactic and extragalactic neutrinos with energies higher than a few 10 GeV. The detector 
consists of 80 strings, each carrying 60 photo- multiplier tubes (PMT) which are being de- 
ployed between 1450 m and 2450 m depth. In addition, there are six more strings with 
a denser array of PMTs between 1760 m and 2450 m, which with the surrounding strings 
constitute the DeepCore detector at the far depth of the IceCube. Since the PMTs at the 
DeepCore are closer, the sensitivity of DeepCore is higher than the rest of IceCube and its 
detection threshold is lower. As explained later, the background events at the DeepCore are 
also highly suppressed such that, in spite of the small volume of DeepCore which reduces 
the statistics of the signal, the discovery chance can still be comparable to that by the whole 
IceCube. 

For neutrinos with energy 0(100 GeV), IceCube cannot distinguish between the neu- 
trino flavors completely. Practically, IceCube can identify two types of events: /i-track and 
shower-like events. Let us explain them one by one. The /x-track events are the Cherenkov 
radiation collected by the PMTs when a muon passes through the instrumented volume of 
the detector. The muons can be produced in two ways: i) from the CC interaction of and 
v^] ii) from CC interactions of Vr (and v-r) which produces the tau (and the anti-tau) particle 
and the subsequent t~ — )• ^~i>^Vr (and t+ — t- /i+f^jZ?,-). However, the small branching ratio 
of the tau particle's leptonic decay, Br(T — t- ^vu) = 17%, makes the contribution from Ur and 
Vr subdominant. Thus, the number of /i-track events is practically given by the fluxes of 
and incident on the IceCube. Considering the position of the muon production point, the 
//-track events for the whole IceCube can be divided into two categories: i) Through-going 
/i-track events for which the muons are produced inside or in the vicinity of the instru- 
mented volume of the detector; ii) Contained /u-track events for which the production point 




(3.1) 
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is within the geometrical volume of the detector. Similar consideration also holds for the 
muon neutrinos detected at the DeepCore. In fact, tracks originated inside the IceCube (but 
outside DeepCore) and passing through the DeepCore will count as "through-going" events 
for DeepCore and as "contained" events for IceCube. 

The second type of events that can be identified at the IceCube is the shower-like events 
which are nearly spherical volumes of Cherenkov radiation from the hadronic or electromag- 
netic cascades. The interactions contributing to the shower-like events consist of: i) Neutral 
Current (NC) interaction of all the three neutrinos flavors z/q and P^; ii) CC interaction of 
V(, and z/g; iii) CC interaction of Ur and i>r, and the accompanying hadronic decay of the 
produced tau particle. The event rate calculation of these interactions can be found in [15]. 

4.2 Observable quantities at neutrino telescopes 

In this subsection, we quantify the observable quantities that we shall use in our subsequent 
analysis; namely, the seasonal variation and the ratio of /i-track events to shower-like events. 
We define two versions of the observable A(ti, Ati; t2, At2) in eq. (2.1). 

• The A^'^(ti, Aii; 5 ^^2) is the seasonal variation of the through-going /x-track events 
for the whole IceCube. To compute A^*-", we calculate the number of ^u-track events 
(per unit of time) with the following formula 



are the scattering CC partial cross section of muon neutrinos off proton and neutron 
respectively, pp ~ 5/9 A^a cm~^ and p„ ~ ^^/S^Na cm~^ are the number densities of 
protons and neutrons in the vicinity of the IceCube in terms of Avogadro's number 
Na, RfM is the muon range in the ice [16], and d = 1 km is the size of IceCube. ^e//) 
which is given in appendix of [17], is the IceCube's effective area of muon detection 
and depends on the zenith angle of the arriving neutrinos 0[t]. Finally, The Sthr is 
the energy threshold of muon detection at the IceCube which we have set equal to 40 
GeV [17]. Notice that in writing eq. (4.1), we have taken into account the fact that the 
maximum equals Ei^^. 

• The second version of the observable defined in eq. (2.1) is ADC(t^,Aii;t2,At2) cor- 
responding to the seasonal variation of the DeepCore /x-track events. The number of 
//-track events determining A^^ is given by 







''f){E^)dE^dE,^ + {u^^v^) 



(4.2) 
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where is the effective volume of the DeepCore given by [5] 

Vj}^{Ef,) = (0.0056 log £^ + 0.0146) 9(275 - E^) + 0.0283 e(£;^ - 275), (4.3) 

in which 0(2;) is the Heaviside step function and E^ is in GeV unit. We have set also 
the energy threshold of the muon detection at the DeepCore equal to 40 GeV [18]. 

In order to employ the shower-like events in our analysis, we define the following ob- 
servable 

^ Number of u-track events , . 

R = — . (4-4) 

Number of shower-like events 

Similar to the observable A, we define two quantities R}*^ and R^^ corresponding to the 
through-going and DeepCore events respectively. In the case of i?^*^, the numerator of eq. (4.4) 
will be calculated using eq. (4.1) and for the denominator we use the geometrical volume of 
the whole IceCube detector. Similarly, for the calculation of the numerator of i?^^ we 
use eq. (4.2); while in the denominator we use the geometrical volume of the DeepCore. 
Generally, if E^^°'"'^'^ > rriDM, where is the energy threshold for the detection of 

shower-like events, there will be no shower-like events at the detector. To carry out this 
analysis, we set equal to 100 GeV. This assumption is rather optimistic as in reality 

^ttr''' ~ 1 TeV. However, in this paper our main focus is on A which does not employ 

^thr 

motivation for pursuing the measurement of A. 



shower-like events. If no practical way to improve 

^shower found, there will be even more 



4.3 Background 

In the measurement of both observables A and R at the IceCube, the background events 
should be considered. The two main sources of backgrounds at the IceCube are atmospheric 
muons and atmospheric neutrinos. The atmospheric muons are isotropically produced in the 
collision of the cosmic rays with the atmosphere. However, from March equinox to September 
equinox (corresponding to spring and summer in the northern hemisphere) when the Sun is 
below the horizon, the Earth acts as a filter for the atmospheric muons so these backgrounds 
will be suppressed. Thus, for the through-going events at the IceCube we integrate eq. (4.1) 
over the March to September period of the year. For the DeepCore, the instrumented volume 
of the IceCube surrounding it acts as a veto (up to one part in 10^ [18, 19]) so it is possible 
to take data at DeepCore throughout the whole year. Thus, we calculate the number of 
muon-track events at the DeepCore by integrating eq. (4.2) over the whole year. The other 
source of backgrounds, atmospheric neutrinos, can be suppressed thanks to the high angular 
resolution of the IceCube in the /x-track reconstruction (about 1° [19]). Considering only the 
events in a cone with half angle 1° around the position of the Sun will reduce the number 
of through-going atmospheric neutrino backgrounds to ~ 6 events per year ( ~ 3 yr~^ for 
DeepCore). There is also an irreducible background from the solar atmospheric neutrinos 
but the fiux is expected to be low [3]. 

5 Numerical calculations 

As mentioned in section 1, the neutrinos at the center of Sun are monochromatic with energy 
equal to the DM particle mass. However, due to the NC and CC interaction of the neutrinos 
with the Sun medium, the spectrum of the neutrinos emerging from the surface of the Sun, 
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in addition to the sharp hne of the monochromatic neutrinos, will contain a continuous tail 
corresponding to the scattered or regenerated neutrinos. We have written a Mathematica 
code to numerically compute the flux of neutrinos at the surface of the Sun. In the code 
we have taken into account the oscillation of the neutrino flavors in the Sun medium, CC 
and NC interactions of the neutrinos with the nuclei, and the neutrinos coming from the tau 
regeneration effect. 

Flavor oscillation of the neutrinos in the Sun medium has been computed using the 
updated matter density profile of the Sun ne{r) from ref. [20]. The oscillation Hamiltonian 
is ^ 

-Hose = ^^-^ + diag{V2GFne, 0, 0), (5.1) 
where is the neutrino energy, Gi? is the Fermi constant and 

(M2)diag = U^M^U = diag(-Am22,0,eAmi2). 

Ml, is the neutrino mass matrix and U is the PMNS unitary mixing matrix. The parameter ^ 
takes values +1 and -1, respectively corresponding to the normal and inverted neutrino mass 
schemes. For the mixing parameters we insert Amf2 = 8.0x 10~^ eV^, Am|2 = 2.4x 10~^ eV^, 
di2 = 33.46° and 6*23 = 45°, 6*13 = and 7°, 6 = 0, tt/2 and tt. 

For the NC and CC interaction cross sections we have used the tabulated values from 
ref. [21]. The NC interaction of a neutrino flavor lowers the energy of the neutrino without 

changing the flavor. The CC interaction of ^I'e and i'^^ transforms them into the corresponding 
charged leptons and therefore reduces the neutrino flux. The tau and anti-tau particles 
created by the CC interactions of Ur and z/^ promptly decay and reproduce neutrinos ("tau 
regeneration"). For example, the decay of the tau particle through the channel r~ — )• ^"v^Vt 
reinject and v-^ with lower energy into the flux of neutrinos. We consider the following 
decay channels for the tau particle: Br(r — eUei'T) = 0.18, Br(r — fiv^v-j-) = 0.18, Br(r — 
TTv^) = 0.12, Br(T aiUr) = 0.13, Br(T pi^r) = 0.26 and Br(r u^X) = 0.13. As can 
be seen, tau regeneration couples the evolution equations of different flavors together (it also 
couples the neutrino flux evolution to the anti- neutrino evolution). The details of the tau 
regeneration calculation as well as the neutrino spectra in the above decay channels can be 
found in [22]. The set of coupled evolution equations describing the neutrino propagation 
can be found in [9]. 

Let us take the state of the neutrinos at the production point to be Ifo). By solving 
the evolution equations numerically, we obtain the state of the neutrinos | Va ] surface) at 
the surface of the Sun. After leaving the Sun, neutrinos will propagate through the empty 
space towards the Earth. The evolution of the state | Va surface) from the Sun surface to 
the Earth is given simply by multiplying each mass eigenstate with phases, exp[imfL/E,y]. 
These phases depend on the distance between the Sun and Earth (i.e., L) which is a function 
of time. Considering the eccentricity of the Earth's orbit around the Sun and also the tilt 
of the Earth's rotation axis with respect to the ecliptic plane, we have calculated 6{t) and 
L{t) during a year (see eq. 2.1 and 2.2). After calculating the values of 6{t), L{t) and the 
spectrum of neutrinos arriving at Earth during the year, we have computed the observables 
R and A according to the formulation presented in section 4. 

It should be mentioned that the flux of neutrinos arriving in the Earth has been cal- 
culated previously by Blennow et al. [8] for different annihilation modes of DM and the 
mixing parameters ^13 = 0, 5°, 10° and (5 = 0. However, since we needed the flux of neutrinos 
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for other values of the mixing parameters (especiahy 5 = 0, 7r/2 and vr), we have calculated 
it independently. With the same input values of the mixing parameters, our results are in 
complete agreement with the results of [8]. 

6 Analysis of information from seasonal variation 

Let us take the amplitude of DM + DM —> 2/0,+ ^ to be Ma/s- The state resulting from the 
DM pair annihilation will be a "pure" coherent state of type \ip) = ^ai3 ■^cii3\'^a{pi) ^(^2))- 

The density matrix is = J^ap^ya Pal3,-yaWa ^^\){T^-t I in which Pai3,^cj = MaijM*^. 

Since the DM pairs are almost at rest inside the Sun, the momenta of the produced pair will 
be opposite to each other (i.e., pi ~ — P2) which means at most one of the emitted neutrinos 
will reach the earth. As a result, the density matrix of the neutrinos pointing towards us will 
be the following "reduced" matrix 

Pap\l'a){l^l3\ with Paj3 = ^ /Oa7,/37 = (A^A4^)a/3 • 

7 

Notice that although the two particle state \tjj) is pure (i.e., plogp = 0), the neutrino states 
reaching the detector are not in general pure (i.e., plogp 7^ 0). In case that the annihilation 
is lepton number violating DM + DM — )• z^^ + z/^, A^^^ is obviously symmetric. As a result, 
p yields Map, i-e-, for p = V pdiagV'^ where F is a unitary matrix, M. = V{p(iiagY^'^y'^ ■ In 
case that the DM annihilation is lepton number conserving (i.e., DM + DM — )• Va'i^p), up 
to subdominant (or perhaps zero) CP-violating effects, we can write Map = -^pa- Under 
this condition, Aia/s can be again derived from p: for a Hermitian Ad, p = YpdiagV^ yields 
M. = V{pdiagY^'^y'^ ■ From Map, information on flavor structure of the couplings of DM can 
be derived. By studying the neutrino flux there is however no way to figure out whether the 
DM annihilation is lepton number violating or conserving. In the following, we discuss the 
possibility of constraining p from the observable quantities like R and A. 
The density matrix p can be diagonalized as 

P = ^na\v'a){v'a\ (6.1) 

a 

where a = 1, 2, 3, > and {i^'aWp) = ^ap- Notice that v'^ can in general be different from 
mass eigenstates z^i, V2 and v^, a-s well as from flavor eigenstates T^fi and Vr- 

Once arrives at the surface of the Sun, it will be a different state that can in general 
be written in terms of the vacuum mass eigenstates as 

surface) = aai|l) +aa2|2) -|-ao3|3) . (6.2) 

We do not a priori know what is the flavor composition of u'^, so aai are unknown. A general 
neutrino state \tp) produced at the Sun center, while crossing the Sun, evolves into IV'; surface) 
as follows 

ts 

1-0; surface) = VF(tp, is) IV) with W{tp,ts) = J\T[1 + iH{t)dt\ 

tp 
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where tp and ts are respectively the production time and the time that neutrino reaches the 
surface and T denotes time ordering. Regardless of the time dependence of H, the evolution 
matrix is unitary so for any arbitrary states {tp) and \x) 

(x; surface! ^/^; surface) = {x\fp)- 

As a result, the states after evolving in time remain perpendicular: 

(f^; surface! z/^; surface) = aajO^j = 5aj3- 

i 

Thus, the knowledge of two rows of the 3x3 matrix, aai, is enough to determine the matrix 

On the way to the Earth, the state will evolve into^ 

!i/;,; detector) = a„i!l) + a„2e*^i^!2) + a„3e*^^^!3) , (6.3) 

in which Ajj = ^mf-L/{2E). Obviously, the number of /x-tracks and therefore both A and 
R depend on X^Q,naP(i^of — ^ ^^J)■ Let us then discuss these oscillation probabilities 

P{T^a ^ I^m) = (6-4) 

i 

As shown in ref. [1], the first oscillatory term given by exp[iAi2] does not vanish as 
we integrate on time over a year when the Earth orbits around the Sun. This term results 
in sizeable seasonal variation. Changing the DM mass, which coincides with the energy of 
neutrino line, this term rapidly oscillates with period 

5m DM _ 47rmgA/ _ mpM ,„ r.\ 

ruDM ~ Aml^L ~ 200 GeV ' ^ ' 

It is noteworthy that even the effects of the phase A13 do not completely vanish. The reason 
is that due to angular momentum conservation, the Earth slows down when it is at the 
aphelion so the eccentricity of the orbit gives rise to non- vanishing average. To quantify this 
claim, let us define 

and 

oisit, At) ^ ^ , ":'':^7,.y"^ - (6-7) 



oi2{t, At) = ^^^^^ . (6-6) 



j;+^\^ArmA,fjit)L-\t)dt 

!l^'''A,ff{t)L-^{t)dt 



For ~ 100 GeV, !0i2! ~ 1 and !0i3! ~ 0.1 which means even the effects driven by A13 do 
not completely average out. The average oscillation probability can be written as 



Jt ^effL \t) 



^The Earth matter effect is irrelevant because for this energy range, the effective mixing is negligible [1]. 
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\aai\^\U^i\^ + 2R[al^aa2U^iU;^Oi2] + 0{aliU%On). (6.8) 

i 

Let us now discuss the antineutrinos. The density matrix of antineutrinos at production 
can be written as 

Pal3\l'a){l'l3\ with = ^ P^q,^/? = (7W^A^*)a/3 • (6.9) 

7 

The density matrix p can be diagonahzed as fohows 

a 

Notice that for the lepton number and CP conserving annihilation (DM + DM — t- uu) with 
M = M\ we find p = p* . On the other hand, for lepton number violating annihilation 
(DM + DM —7- V + v,!? + v) with M = A4'^, we find p = p- In any case in eq. (6.1) for 
neutrinos and the one in eq. (6.10) for antineutrinos are the same. Moreover, {v'^) and at 
the source are the charged conjugate of each other. However, the matter effects on and 
1 1^4) are different so surface) will not in general be the charged conjugate of | ; surface) . 
Let us expand the antineutrino state at the surface as 

surface) = Oaill) + aa2\^) + OasIS). (6-11) 

As explained above because of the matter effects on the way from the Sun center to its 
surface, o* ^ is not in general equal to Uai- It is straightforward to confirm that the oscillation 
probability for anti-neutrinos, P{ua — ^ z^/^)) and its average are given by formulas respectively 
similar to Eqs. (6.4,6.8) replacing aai with Uai and U^i with U*^. 

Since the scattering cross sections of antineutrinos are smaller than that of neutrinos 
by a factor of about 2, their contribution to the events at the detector will be smaller. 

Obviously, replacing p with r, we obtain the formula for {P{'Ua — ^ ^r))\\^^^ and 
{P{i>a —5- i>r))\\^^^ . Notice that to the leading order in ^13 and ^23 — 7r/4, \U^i\ ~ \Uri\ 
and UriU*2 — U^iU*2- a result, at this limit, P{i'a — ^ i^t) does not carry new information 
relative to P{i'a t^^i)- Muon-track events dominantly receive contributions from the CC 
interactions of and z^^. A subdominant contribution also comes from the CC interaction of 
Ur and the subsequent decay of tau to the muon which is suppressed by Br(r — )• pi^v) ~ 17%. 
As a result, up to a correction suppressed by Br(T — t- pw') sin^ia, muon track events within 
interval (t, t + At) from the sharp line is proportional to 

K{t,At) = 5]n, ({P{i.^ ^ z.^))|*+^* + ^AP{i^a ^ ^m))!^''*) , (6.12) 

which should be added to the contribution from the continuous part of the spectrum. In 
general, the //-track events can be written as 

A + BK{t,At), (6.13) 

where A comes from the continuous part of the spectrum. A priori, both A and B are 
unknown. While 

(j) = Am\2L / {2mDM) ^ 27r , 
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Figure 1. Number of the through-going /^-track events from the spring equinox (20 March) to autumn 
equinox (22 September) for ^13 = (left) and ^13 = 7° (right). In drawing these histograms, we have 
assumed tliat DM particles with mass tudm — 270 GeV decay via DM+DM— ^ i^e^e- Foi" the DM 
capture rate in the Sun, we have taken Cq = 3.4 x 10^^ s^^. 



number of muon (racks 



rrunrbcr of rrruorr (racks 



Tirrrc of year 

20 Apr 20 May 20 June 20 July 20 Aug 22 Sep 



20 Apr 20 May 20 Jui 



Tinre of year 
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Figure 2. The same as figure 1 except that DM+DM— > v^v^ 



the variation in (j) due to seasonal change in the Sun Earth distance AL ~ 5 x 10^ km is 
of the order of 27r. This means by studying the seasonal changes in muon track events, the 
ratio Am\2/fnDM can be derived. 

Histograms in figures 1 and 2 display such seasonal changes for the case that DM pair 
dominantly annihilates into neutrino pairs. The horizontal axis shows the months between 
spring and autumn equinoxes during which neutrinos pass the Earth to reach the IceCube. 
As discussed in the previous section, during these months atmospheric muon background 
due to absorption by Earth is reduced and the data is therefore reliable. The peaks and dips 
in the histograms of figures 1 and 2 respectively correspond to the peaks and dips of the 
oscillation probabilities Pive i^^i) and -P(f^ — )• i'^). If the present bound on the neutrino 
flux is saturated, after ten years of collecting data, the number of events during each month 
can amount to a few hundred so extracting Am^g/'^^DA/ with a reasonable accuracy (about 
10 %) would be possible. Since Amf2 is already known, this ratio yields the DM mass. The 
value of the DM mass can in principle be derived by accelerators such as the LHC. If the 
two values coincide, it will be a noteworthy confirmation of the validity of the approach. In 
addition to Am'^ /itidm, the seasonal variation also yields 

A + B^na{\aai\^\U^^\^ + ^\aa^\^\Uf,,\^), (6.14) 
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as well as 



(6.15) 



Moreover, the seasonal variation yields the following combination of 0, arg[a* ]^aQ,2] and 
arg[aaiaa2]: 



arg 



+ </). (6.16) 



Since (/)/(27r) ~ 25(200 GeV/m£)Af), to derive meaningful information on arg[a*-|^ao2] and 
arg[a*|aa2]) (p has to be known with percent level accuracy or better which does not seem to 
be achievable in foreseeable future. To derive these four independent pieces of information 
(i.e., Am^2/"^-DM and three combinations in eqs. (6.14,6.15,6.16)), one can divide the time 
period between two equinoxes to four segments and measure the number of /_f-tracks in 
each interval. If the number of events at each interval exceeds a few 100, the statics will be 
enough to perform the analysis and derive the aforementioned combinations with about 10 % 
accuracy. 

Notice that the density matrix in eq. (6.1) (or equivalently matrix p in eq. (6.10)) 
contains seven free parameters: three real values of ria and four parameters (3 angles + 
1 phase) determining u'^ in the flavor basis. Of course, the three combinations that can 
be extracted are not enough to reconstruct p (and therefore Ma/^)- However, these three 
combinations yield valuable insight on the flavor structure of pap and therefore In 
the following, we discuss the additional information that can be derived from the shower-like 
events. 

As is well-known, shower-like events receive contribution from three sources: (1) uni- 
versal NC interaction of all neutrino species which is insensitive to neutrino oscillation prob- 
ability because Yl,pP{^oi — ^ ^p) = Y^p^i^a — ^ 1^/3) = 1; (2) CC interactions of z/g and 
whose contributions are given by P{i'a — ^ J^e) + Pi^^a — ^ i^e)f(^'e)/c(fe); (3) CC in- 
teractions of Vr and v'r and the subsequent hadronic decay of r whose effect is given by 
[1 — Br(r — )• pi'i')][P{ua — )■ z^r) + P{i^a — ^ i^r)'7(^'T)/<7(^'r)]- Remember that up to a cor- 
rection of the order of Max[0i3,023 — P ^t^I^^t^^) — -P(W— ''W\ One can write the 
shower-like events as 

C -rl3K{t,At) , (6.17) 

where K{t,At) is defined in eq. (6.12) and r is a calculable quantity (see section 4). The 
quantity C receives contribution both from continuous part of the spectrum and the sharp 
line. Combining the number of shower-like events with the information derived from the 
seasonal variation of the /x-track events yields C + rA. In principle, C can be predicted 
from A and B but for such a prediction, knowledge of the shape of the spectrum is required 
which to some extent can be extracted from energy measurements. In fact, for a given 
value of m£)M, the shape of continuous spectrum from various annihilation modes such as 
DM -|- DM —7- bb,Tf, ZZ,W'^W~ can be determined. The shape of spectrum for scattered 
and regenerated neutrinos from DM -|- DM — t- z^i/ can be also determined. In other words, 
C depends on the relative ratios of the different annihilation modes. Such an analysis might 
provide a cross-check for the consistency of the analysis and the measurements. For further 
study of the role of R in determining the decay modes see ref. [24]. 
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7 Illustrative examples 



From the theoretical point of view, the neutrino states produced inside the Sun do not need 
to coincide with any of flavor eigenstates as the flavor structure of the physics governing the 
DM sector might be different from the physics determining the flavor structure of the SM. 
However, for illustrative purposes in this section, we study the case that the neutrino state 
at production corresponds to u^, or u^. In other words, we assume p and p at production 
to be flavor diagonal. Remember that the study applies both to DM + DM — )• fa + i^a^T^a + ^a 
and DM + DM z^q, + Va- 

At the production point, the matter effect \plGpne is much larger than Am^2/(2-E') 
so It'e) and |z?e) will respectively correspond to \v2] and \i>\). As propagates outside, it 
will reach the 12-resonance region. At these energies, the transition is non-adiabatic so at 
the surface |aei| ~ |'2e2|- For the case of antineutrinos, there is no such resonance region 
so at the surface |aei| ^ |ae2|- For normal hierarchical scheme and nonzero ^13, neutrinos 
will undergo a second resonance. For ^13 within the reach of the forthcoming experiments 
[23] (i.e., ^13 > 7°), this resonance is also non-adiabatic. For inverted hierarchical scheme, 
instead of neutrinos, antineutrinos will undergo the resonance due to the 13-splitting. 

For lower values of the DM mass, detection rate becomes smaller as the neutrino-nuclei 
cross section decreases with lowering the energies of neutrinos. On the other hand, with 
increasing tudm, the sharp line composed of the unscattered neutrinos diminishes, leading to 
a suppressed A. It can be shown that for ^ 500 GeV, the tail of the spectrum from the 

scattered neutrinos suppresses A down to values A < 0.5. Thus, generally we expect that 
in the range, 100 GeV < ttidm < 500 GeV, the A measurement method will be most useful. 
We focus on m/) ~ 200 GeV and study the high sensitivity of A to moM. We have checked 
for other values of rriDM in the range 100 GeV to 300 GeV and found that our results are 
robust against varying moM- Using the code described in section 5, we have calculated the 
seasonal variation, Rdc and iZthr for niDM = 197, 200, 203 GeV, ^13 = 0, 7° and 5 = 0, 7r/2, vr, 
taking the initial neutrino flavor to be v^, and v^- The results are displayed in tables 1, 2 
and 3. 

For the shower- like events, the effective volume of a detector is approximately equal to 
its geometrical volume: Vg^j"*''' = dA^ff where d and A^jf are respectively the depth and 
effective area of the detector. However, since the muons can penetrate farther, the effective 
volume for their detection can be larger: V^jj = {R^ + d)Af.ff where i?^ is the muon range 
which depends only on the chemical composition of the medium. As a result, the ratio of 
//-track to shower-like events is enhanced by a factor of (R^ + d) /d. Obviously, for a smaller 
detector this enhancement factor is larger. That is why the ratio in the case of Deepcore, 
Rbc, is so much larger than the ratio in the case of the whole detector, i?thr- (See tables 1, 
2 and 3.) 

The matter effects inside the Sun do not distinguish between and I'r- Thus, if the 
neutrino mass matrix is /ir-symmetric in vacuum (i.e., ^13 = 623 — 11/4: = 0), the neutrino 
sector remains /ir-symmetric in the matter. This means that in this limit, regardless of 
whether neutrino mass scheme is normal or inverted, we can write 
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Table 1. Seasonal variation and //-track to shower- like ratio for DM particles with mass m£)j\/ (in 
GeV) and annihilation mode DM -|- DM v^v^. N/I indicates normal versus inverted neutrino mass 
scheme, i?^^ indicates the ratio of the numbers of muon-track to shower-like event for through-going 
events measured by the whole IceCube. BP^ is the same quantity measured by DeepCore. A20Mar 
is the seasonal variation between two equinoxes A20Mar = A(20 Mar, 186 days, 23 Sep, 179 days). 
Finally, AgApr = A(3 Apr, 186 days, 6 Oct, 179 days). See eqs. (2.1,4.2,4.4). 



and 

|z.,;surface) = ^^|1) - ^^|2) + ^|3) . (7.3) 

^/2 

Notice that is different from Qyi mixing angle in the ?7mns neutrino mixing matrix. 
Similar relations holds for antineutrinos but with a different value of Q\2- For antineutrinos, 
we expect cos^f2 — 1. In terms of the notation in eqs. (6.2,6.11), = cos0f2, aei = cos(9f2, 
ffleS = Oe3 = 0, 0^3 = 0^3 = 0^3 = Ourz = l/\/2 and etc. From the above relations we find 

I aji I = I Or! I, a*^iaf,2U^iU*2 = a*^iar2U t,iU*2 and a*^ia^2U*iUfj,2 = a*i^T2U*iU^2- Notice 
however that a*^iafj.2Ufj,iU*2 / a*iaT3U^iU*^. Thus, in the /ir-symmetric case, up to the 
corrections of order of O13 (see eq. (6.7) for the definition), (-P(f^ — )■ t'^)) ~ {P{i^t — ^ i^^i)) and 
{P{v^ —7- v^)) ~ {P{vt — )• v^)). This means, to leading approximation, if we take only the un- 
scattered neutrinos, the contribution to A and R will be the same for initial and u-r. This 
makes distinguishing between (DM+DM — >• I'^i'^) and (DM-I-DM — t- VtT^t) difficult. However, 
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Table 2. The same as Table 1 for DM+DM^ v^,v^,. 



the regenerated neutrinos coming from i/,- — ?■ r — t- i^t- in the Sun differentiate between them. 
Remember that at these energies Aml^RQ/2E < 2tt. In fact, fn can only partially oscillate to 
I'r before exiting the Sun. Thus, the contribution to regenerated neutrinos strongly depends 
on the initial flavor of the neutrino. Of course the regenerated neutrinos will have less 
energy, so they are more important for DeepCore than the whole IceCube because the energy 
threshold of DeepCore is lower. Our results show that for ^13 = 623 — 7r/4 = 0, the deviations 
of A^*-" and E}^ for initial from the same quantities for initial Ur are negligible (compare 
tables 2 and 3). As seen from the tables, this deviation for the case measured by the DeepCore 
(i.e., A^'~' and i?^*^) is larger. This can be explained by the contribution from the regenerated 
neutrinos. 

Another point to note is that since the validity of eq. (7.1) does not depend on scheme, 
we expect that for ^13 = and initial fgi A and R to be the same for normal and inverted 
mass schemes. Results shown in table 1 confirm this expectation. 

The tables show that the sensitivity of A to mi^M is very high. As discussed in the 
previous section and shown in figure 1, by studying the variation over four intervals of time, 
Ami2/mDM can be derived. 

From the tables we observe that for (DM+DM — )• Ve+Ue), A can be significant; however, 
for the cases (DM + DM i^^ + f^) and (DM + DM z/^ + z/^), A is in general smaU. We 
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Table 3. The same as Table 1 for DM+DM^ v^Vr- 



found that this observation is robust against varying m£)M- The reason is two folded: (i) 
For initial and v-r, the contribution from the continuous regenerated neutrinos to the flux 
at the detector is higher so the denominator of the ratio in eq. (2.1) is enhanced leading 
to lower values of A. From tables 2 and 3, we observe that A^^ is smaller than A^^. The 
reason is that DeepCore has a lower detection threshold so it receives a larger contribution 
from regenerated neutrinos which are in general less energetic, (ii) The second reason is that 
for a = /i, r; aaiaa2U^iU^2 is smaller (see eq. (6.8)). However, we should remember that 
(DM + DM —7- f^ + f^) and (DM + DM — )• Ut + i't) are quite specific cases. In general DM pair 
can annihilate to general coherent combinations of the different neutrino flavors for which A 
is in general large. 

8 Conclusion and discussion 

We have shown that by studying the time variation of /i-track events, valuable information 
on the nature of DM particles can be derived: (i) Measuring a nonzero value for A defined in 
eq. (2.1) implies that there is a sharp monochromatic feature in the neutrino spectrum which 
in turn means DM + DM v + u is [one of] the significant annihilation modes, (ii) Once 
such a feature is established, by studying the /x-track events over four different time intervals. 
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the value of Am^g/mDM as well as three independent combinations (see eqs. (6.14,6.15,6.16) 
as well as eq. (6.13)) of aai and (iaj which are defined in eqs. (6.2,6.11) can be extracted. 
Thus, by this method, moM can be extracted and checked against the value derived from 
studying the endpoint of the spectrum or from the accelerator experiments. Although the 
information will not be enough to completely reconstruct the flavor structure of the amplitude 
of DM + DM — 7- f Q. + Vji, the combinations in eqs. (6.14) and (6.15) contain invaluable 
information that can constrain models predicting DM + DM — )■ t'a + Vfj. In particular, 
observing nonzero oscillatory effects rules out models predicting democratic neutrino flavor 
production. Moreover, our numerical results show that by studying A, models predicting 

DM + DM —7- W can be distinguished from those predicting DM + DM — t- or 

DM + DM —7- Ut-+ ^u} . Notice that to derive information on AAap from aai, the evolution of 
the neutrino state from the Sun center to its surface has to be taken into account. This can 
be readily done by codes such as the one developed to derive the numerical results of the 
present paper (see, tables 1, 2, 3). 

We showed that it is possible that while DM particles dominantly annihilate into neu- 
trino pairs, their interaction rate with nuclei in the Sun is high enough to lead to high 
enough DM capture rate and therefore to a neutrino flux saturating the present bounds. 
If the present bound on the neutrino flux from dark matter annihilation is saturated, after 
about 10 years of data taking, the statistics will be high enough to make this method reliable. 
In this method, only /i-track events are employed so high energy threshold for the detection 
of shower-like events will not be a hinderance. To carry out this analysis, it will be enough 
to take data from 20th of March to 23rd of September when the neutrino from the Sun pass 
through the Earth to reach the detector and as a result, the atmospheric muon background 
will not be a problem. The method presented in this paper can be combined with the energy 
spectrum measurements [6] and the ratio of /i-track to shower-like events [24] to draw a more 
complete picture of the dark matter annihilation modes. 
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